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Abstract: We study nuclear symmetry energy of dense matter using holographic QCD. 
We calculate it in a various holographic QCD models and show that the scaling index of 
the symmetry energy in dense medium is almost invariant under the smooth deformation 
of the metric as well as the embedding shape of the probe brane. We find that the scaling 
index depends only on the dimensionality of the branes and space-time. Therefore the 
scaling index of the symmetry energy characterizes the universality classes of holographic 
QCD models. We suggest that the scaling index might be also related to the non-fermi 
liquid behavior of the interacting nucleons. 
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1. Introduction 

Nuclear symmetry energy is one of key words in nuclear physics as well as in astrophysics. 
Its density dependence is a core quantity of asymmetric nuclear matter which has important 
effects on heavy nuclei and is essential to understand neutron star properties. A big surprise 
is that such an important quantity is still poorly understood after 80 years of its definition, 
especially in the supra-saturation density regime. See references [2|, [|, ||, ||, |(| [j], ||. |9| 
for a review and for a recent discussion. Not much data is available from experimental 
side and theoretical calculations showed all possible results in high density such that no 
consensus could be reached: some showed stiff dependence on density, while others showed 
soft one depending on models and parameters. 

Given this situation, it would be very interesting if we can examine the behavior of 
the nuclear symmetry energy at high densities with a reliable calculational tool. Recently 
|1C] we used a gauge/gravity duality jO], 12, [H| |l4], 15] to calculate the nuclear symmetry 



energy. We treated the dense matter in confined phase using the method developed in our 
previous paper fllq , |l7]] . Our result showed that the symmetry energy should be stiff in high 
density and its density dependence goes like ~ p 1 / 2 . There, we attributed the stiffness to 
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the repulsion due to the Pauli principle and suggested the relation of the scaling exponent 
to the anomalous dispersion relation. 

The purpose of this paper is to examine how universal or robust is the result. If the 
result changes under small variation of the gluon dynamics, the result is not very interesting 
since the true QCD dual is not yet known. Only when the result is largely background 
independent, it can be considered as an interesting one. The universality of the rj/s is the 
reason why it is interesting even though it is not calculated in the QCD itself. 

We will first calculate the deformation of the metric under certain class of the D brane 
configurations and show that the result is not much dependent on the metric deformation. 
The scaling behavior is rather insensitive whether we use the flat embedding or exact shape 
of the brane embedding, showing the universality of the result. On the other hand, we will 
see that the scaling exponent depends on the dimensionality of the color and flavor branes. 
We call such discrete dependence of the scaling dimension as the universality class of the 
symmetry energy. 

The rest of the paper is summarized as follows. In section 2, we give a definition 
of symmetry energy and general formula in the brane set up. In section 3, we calculate 
the symmetry energy S2 in nuclear matter using D4- as well as D3- confining geometry 
for various probe branes. In section 4, we calculate £2 in quark matter using D4, D3 
de-confining geometry. In section 5, we reproduce the scaling exponent analytically using 
the BPS background and flat embedding and thereby argue that it is a invariant under the 
smooth deformation of the gluon dynamics. In section 6, we discuss the possible relation 
of the scaling exponent with the non-fermi liquid nature of the strongly interacting dense 
matter system and conclude. 



2. Symmetry energy 

Let N, Z be the nucleon and proton numbers respectively, A be their sum and p is total 
baryon density. Then energy per nucleon in nuclear matter system can be expanded as a 
function of the isospin asymmetry parameter a = (N — Z)/A, 

E(p, a) = E(p, 0) + E sym (p)a 2 + 0(a 4 ). (2.1) 

The bulk nuclear symmetry is defined as the coefficient E sym (p) in the above expansion. 
There is no term which is odd power in a due to the exchange symmetry between protons 
and neutron in nuclear matter. It is a energy cost per nucleon to deviate the line Z = N. 

To calculate nuclear symmetry energy in holographic QCD, we introduce two flavor 
Dq branes for up and down quarks in the metric background created by the -/V c of Dp color 
branes. For simplicity, we assume that masses of up and down quarks are the same so that 
two branes have same asymptotic positions. For the confining geometry, we can introduce 
a baryon by the baryon vertex [18] which is a compact D(8 — p) branes wrapping the S 8 ~ p 



transverse to the D(8 — p). If there are Q of them we distribute them homogeneously along 
the 3 non-compact spatial direction of Dp. From each baryon vertex, N c of the strings 
emanate and end at one of the probe branes. Let Q\, Q2 strings end on up and down 
branes respectively. The end points of the strings have £7(1) charges that will create the 
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U(l) gauge field on each brane. Such U(l) charges are responsible for the quark density of 
each type of quarks. 

The total free energy of the system can be written as 

•Ftotal(Q) = HQ) + ?D q (Ql) + ^D q (Q2), (2-2) 

where Q is number density of source and J~(Q) is a quantity which depends only on total 
charge Q as we will discuss later. We can define total charge density and asymmetry 
parameter as 

Q = Q 1 + Q 2 , a = 9l_9l. (2.3) 

If we fix the asymptotic value of two probe brane to be same, the total free energy has 
minimum at a = 0(l7]]. Then we can expand total free energy in a; 

•Ftotal(Q) = E + E 1 a + E 2 a 2 + ■■■. (2.4) 

The first term, Eq = J-{Q)+2J r D q (2J 1 can De identified with the free energy for symmetric 
matter. The second term in (|2.4| ) is zero because ( |2.2D is symmetric in Q\, Q 2 . The 
symmetry energy is defined from the energy per nucleon and given by 

e 2 (Q)_q a^Qi 



5 2 



Q 4 dQi 



(2.5) 

Q/2 



To calculate symmetry energy from D-brane set up, we consider probe Dq brane spans 
along (t,Xd,p) and wraps S n where n = q — d — 1. The induced metric on Dq brane can 
be written in general form; 



dsjjg = -G tt dt 2 + G xx dxj + G pp dff + Gnndfi,_ d _i. (2.6) 

To introduce number density, we turn on time component of U(l) gauge field on the probe 
brane whose action is given by 

S Dq = fjL q j da q+l e-*^det(g + 2na'F), (2.7) 

where fj, q is tension of Dq brane. The free energy can be identified with the action for fixed 
charge sector, which is the Legendre transformation of the original action with respect to 
the gauge field, 

-WQ) = r q j dpy/G^ p y]& + e-^G&Gftn, (2.8) 

where r q = p q Vd^l n , Q = Q / {^a' T q ) and n = q — d — 1. From ( |2.5|) we symmetry energy 
is given by; 

S 2 = 2r q I dp QV^Gy P e-^G i^_ ^ g) 



For details, see Appendix |A[ This result can be applied for general Dp brane background 
and we will apply it to confined phase as well as deconfined one. 
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3. Symmetry energy of nuclear matter 

In this section, we will discuss symmetry energy in the confined phase. There are several 
examples of the metric background corresponding to the confining phase. In this paper, 
we will consider two examples based on DA and D3 branes. In the case of DA brane, the 
geometry is obtained by (double) Wick rotating time and a compact spatial direction. In 
D3 brane case, we will use the non-supersymmetric geometry with nontrivial dilaton field 
|19| j. In such geometries, the net force on the probe brane is repulsive. 

We identify the the chiral symmetry as the the symmetry rotating the probe brane 
in the transverse plane following j2(|. In the limit of zero quark mass, this symmetry is 
spontaneously broken due to the repulsive nature of the net force. As a consequence, the 
value of chiral condensation has finite value and the baryon vertex is allowed to exist. The 
baryon vertices play a role of source of U (1) gauge field on probe Dq brane as we discussed 
in previous section. 

To make the system stationary, we have to impose 'force balance condition' between 
baryon vertex and probe brane. The details of the solution of baryon vertex and force 
balance condition are in Appendix [B[ The symmetry energy can be understood as the 
energy costs when the system is deviated from symmetric matter. To achieve the deviation, 
we need to attach different number of charges(strings) on each brane, which in turn gives 
different embedding for each probe brane. The symmetry energy (|2.9|) is nothing but the 
energy difference of the D-brane system between symmetric and asymmetric distribution 
of source on probe branes. The schematic figure is drawn in Figure. [l|. We first solve the 
equation of motion for probe brane numerically with given charge and quark mass together 
with force balance condition. Then substituting the solution to ( |2.9| ), we can calculate the 
symmetry energy. 




Figure 1: Schematic figure of D-brane configuration for (a) symmetric matter and (b) asymmetric 
matter. 



3.1 DA brane background 

Here we consider probe D6, DA and D2 branes in the DA background which is given by 
ds 2 m = (|) ^ (r^cteW + f(U)dxl) + (^) 3/2 (J^ + UHSII} 
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£4, /(C/) = l- (^) 3 , ^ = 77^/3. (3.1) 



ft. 



The Kaluza-Klein mass scale Mkk is defined as inverse radius of the X4 direction: Mkk - 

I ^ fl . The bulk parameters, Ukk, 9s, R and the gauge theory parameters Mkk, 9ymi ^ := 
g YM N c are related by 



A 



2nl s N c M KK 



U KK = t:\MkkII R 3 
y 



2M KK 



(3.2) 



We introduce new coordinate £ by -4- 



^4 = (|) + ^ + W)dxl) + ^ (|) ' («£» + f cfflg) • (3.3) 

The relation between U and £ is 

(J3/2 = £3/2 



and obtain, in Euclidean signature, 



3/2 



1 + 



£3/2 



L0+, 



(3.4) 



where JT 3 ^ = 2£ 3 ^. We res cale £ = £,kk£ such that singularity is located at £ = 1. Then 



the background (3.3) and the dilaton can be rewritten as 



( \ 3/2 1/3 

<i«i , ( ^ ) £ 3/ V (dt 2 + df 2 + /(t/)^ 2 ) + (i? 3 a„) 1/2 ^ (a, 2 + £ 2 ^i) , 



R 



(3.5) 



The baryon vertex in this background is a spherical DA brane wrapping S 4 . The induced 
metric on compact DA brane is given by (see |B.2j for the detail), 



ds%y = (^) 3/ ^ 3 /Wt 2 + (i? 3 ^A-) 1/2 -f ^ /2 



The free energy of the compact DA can be written as 



1 + 



4tt I dO 2 + sin 2 ^ 2 



. (3.6) 



with 



•Fsy = r 4 y d0yJu%*(P + £ ,2 )/W + sin 6 



5(0) = -2 + 3(cos - - cos 3 0). 



(3.7) 



(3.8) 



The embeddings of baryon DA brane can be obtained as a function of £0, the position of 
DA brane at £ = 0. From the embedding solution, we can get value of £ and its slope at 
9 = 77 which will be used in force balance condition. The details of the embeddings are 
discussed in 
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Now we add probe brane in this geometry; we can put D6, DA and D2 brane corre- 
sponding to the 3 + 1, 2 + 1 and 1 + 1 theory at the boundary. The general induced metric 



on Dq(q = 6, 4, 2) brane can be written from ([A^ 

3/2 



ds 2 



D a 



jKK 

R 



e /2 ^ + {dt 2 + dx 2 ) 

1/3 

3* W2^+ 



+ (R 3 Ikk) 



£3/2 



l + Y^dp' + p'dn^ 



(3.9) 



From ( A.10| ), we can write the free energy of Dq brane as follows 



?d 9 = r q I dp^uf{l + y 2 WQ 2 + P 2n oJ 



(3.10) 



where n = q — d — 1 and r q = T q ^ K . Q is related to the number of source Q as follows 



Q 



Q 



T 9 (2 7 ra')4V ) 



(3.11) 



The boundary condition is described in the appendix B ( B.14j ). We now can solve the 
equation of motion with force balance condition for given Q. 



For the nuclear symmetry energy (2^9) in large iV c theory, we need to define what is 
the proton and neutron. In the case of N c = 3, proton consists of two up quark and one 
down quark( uud), and neutron is udd. But in generic _/V c , it is not clear what is the quark 
contents of proton and neutron. There are many possibilities for quark configurations of 
proton and neutron but we consider two possibilities in Figure S 



N c + 1 N c - 1 





(a) (b) 
Figure 2: proton and neutron in generic N c . 



• case (a): In this case, the difference of quark number between proton and neutron 
is always 1. To make this configuration be possible, we assume that iV c is odd. From this 
configuration, we can set 

Ql - Q2 = N p - N n , Q 1 + Q 2 = Q = N B N C , (3.12) 

where N p is number of proton and N n is number of neutron. Then, a can be written as 

. Q1-Q2 N p -N n 

a= Q^Q- 2 =^N^' (3 - 13) 
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and 

1 + a 1 — a 

Ql = —^-N C N B , Q 2 = — - — N c Nb- (3.14) 

From this definition, the second order term becomes 

. 2p _ fN p -N n \ 2 „ _ fN p -N n \ 2 E 2 



Then, the symmetry energy per nucleon can be identified as 

From (|2.5|), we can get symmetry energy per nucleon in terms of elements of induced metric, 



Ikk i i- q^W+^) p 2n 4 {d ' 1] 



S2 = ^~ dp^-^- " ■ (3.17) 



There is 1/Nc factor in the symmetry energy ( |3.17 ) which implies that the symmetry 



energy is suppressed by Nc- It is consistent with the definition of proton and neutron: 
there is only one quark difference between proton and neutron and therefore, for large Nc, 
it is not easy to distinguish these two particle and hence symmetry energy becomes zero 
for large Nc- 

• case (b): In this case, proton consist of Nc — 1 up quarks and one down quark, and 
neutron has single up quark and Nc — 1 down quark. The total difference and total number 
can be written in term of proton and neutron number as follows 



Then, 



and 



Qi-Q2 = (N C -2)(N P -N n ), Qi + Q2 = N C N B . (3.18) 



- Q1-Q2 (N c -2)(N p -N n ) 

01 = Q^TQ- 2 = NcW B ' (3 - 19) 



N p -N n \ 2 (N c - 2) 2 
N B J N c 

Performing same procedure with previous case, we get symmetry energy as follows; 



a^ 2= (±>_^) e 2 . (3.20) 



N c 



.. l KK (N C - 2? f Q^f(l + Y*) gcj^ „ 
#2 = — • ~ T / dp— — . (3.21) 

Q 2 + Ap 2n Lof M 



Instead of suppression by Nc, the symmetry energy grows with Nc factor for large Nc 
limit. 

Considering other intermediate case in similar fashion, we can easily see that the free 
energy should be lowest in the case (a). Therefore we take the definition of proton defined 
in (a). 
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We can convert all parameters in terms of physical quantities such as 't Hooft coupling 
A, Kaluza-Klein scale Mkk and density. The nuclear density and quark mass can be 
written as 



Q 



2 (d-l)/3 n 



m„ 



N c V d 32(d-l)( 27r )n+d-l 

2 2 / 3 -9^ ' 



" X d - 2 M d KK Q, 



(3.22) 



and the coefficient in ( |3,17| ) becomes = "^-XMkk- By substituting these numerical 
solution into (2.9) we can get symmetry energy for each embeddings in terms of density 
and quark mass. From the meson mass calculation, we choose 



A = 18, M KK = 1.04GeV, N c = 3. 



(3.23) 



With this values we can calculate density and quark mass dependence of symmetry energy 
for D6, DA and D2 probe brane cases. See Figure y. 




(a) 




D2 brane 



(b) 




(c) 

Figure 3: Density dependence of symmetry energy for several quark mass for (a) D6 probe brane, 
(b) D4 probe brane and (c) D2 probe brane. 



• D6 probe brane: 

In this case, the boundary theory is in 3+1 dimension. From the Figure || (a), the symmetry 
energy with D6 probe brane seems to have square root behavior. All the lines are well fitted 
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to S 2 = S (g/g ) 1/2 , with 27MeV < S < 36MeV. For small quark mass(5MeV), symmetry 
energy curve is fitted to S 2 ~ 28(MeV)(^/^o) 1 ^ 2 - As quark mass increase, the symmetry 
energy curve move downwards, in other words, it become softer up to quark mass is around 
100 MeV. After then, the symmetry energy curve moves to upwards (stiffer) as quark mass 
increases. 

• DA probe brane: 

The dual theory lives in 2+1 dimension and the result is drawn in Figure ^ (b). In this 
figure, most symmetry energy curves are on top of each other unless quark mass is very 
large. And it's behavior seems to be linear for small density. 

• D2 probe brane: 

In this case, dual theory is of 1+1 dimensional. The density dependence in the free energy 



( 3.10 ) factors out: 

T D2 = f 2 J dp^Q 2 + l^Joj 4 J 3 {l + Y 2 ), (3.24) 

so that the embedding configuration is independent of density. The symmetry energy can 
be written as follows 



^KK Q f , I 4/3 



S 2 = ^ „ ^ / dWwr(l + Y 2 ). (3.25) 

4vra' (Q2 + 4)3/2 _/ P V + V > V ) 

The embedding configuration looks flat for any quark mass which means that the symmetry 
energy is almost same for all quark mass region. See Figure. |3| (c). 

3.2 D3 brane background 



We now consider the confining geometry based on D3 brane [21 

ds 2 w = e^ 2 (^A 2 {r)r 1tlu dx^dx u + ^dr 2 + i? 2 ^ , (3.26) 
where A(r) and the dilaton are given by: 

while the five-form remains unaltered from the pure AdS solution. R is the AdS radius 
and ro is the position of the singularity determined by the value of gluon condensation. 
Unlike D4 brane background, the exact relation between ro and < G 2 > is not clear. So 
we will see power behavior of symmetry energy in terms of density and asymptotic value of 
probe brane. For simplicity, we set ro = 1. Procedding similarly with the previous section, 
we can get Hamiltonian for baryon D5 brane; 

F BV = r 5 J d6A(r) ^ e^{r' 2 + r 2 )^J D{6) 2 + sin 8 6 , (3.28) 

where, 

D{6) = -i;0+ 5 sin0 cos 6 + sin 3 6 cos 6. (3.29) 
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The details of this embedding are discussed in p^ ]. We add D7, (or D5, Z)3) brane as a 
probe. The induced metric on Dq brane can be written as 



ds 2 Dq = 



,2 



^A 2 (r)(dt 2 + dx 2 ) + ^ {(1 + Y 2 )dp 2 + p 2 dn 2 q _ d _ 1 } 



(3.30) 



And the free energy for probe brane is 



T Dq = r q j dpyfe<l>A{r) 2 {l + Y 2 )\/q 2 + e^ 5 U( 



j. pin 



(3.31) 



The symmetry energy ( |2.9| ) becomes 



1 f, Qp 2n eM^A(r) 2d+ Wl + Y 2 



Ana' 



dp 



Q2 + 4 Le Ui- 5 )A(r) 2d p 2n ) 



3/2 ' 



(3.32) 



With same method of the previous section, we can calculate symmetry energy, see 
Figure f§. 





(a) 



(b) 




(c) 

Figure 4: Density dependence of symmetry energy for several quark mass. 



• D7 probe brane: 

In this case, the corresponding boundary boundary theory is of 3+1 dimension. The density 
and dependence of symmetry energy is drawn in FigureH. In this figure, one can see that 
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symmetry energy has power behavior unless asymptotic value of probe brane(Yoo = 0.1) is 
very small. Actually, the symmetry energy curves with 1 < < 50 behave as 

4 = ^, (3.33) 

with 0.21 < 5 < 0.31. 

• D5 probe brane: 

The boundary theory is of 2+1 dimension and the symmetry energy behaves 

S 2 - ^ 1/2 , (3-34) 

with 0.34 < So < 0.5 unless Yoo is very small, see Figure § (b). 

• D3 probe brane: 

The boundary theory is of 1+1 dimension. We can see the symmetry energy grows linearly 
for small density, see Figure ||| (c) . 

For the given boundary space-time, the density dependence of symmetry energy seems 
to depend on the dimensionality of D-brane system we use. For the 3+1 dimensional 
boundary theory, the probe brane should be D6 for DA brane background and it should be 
D7 for D3 background. In the case of D6 brane, 5*2 ~ g 1 ^ 2 . On the other hand S2 ~ f? 1 ' 3 
for D7 probe in D3 background. These differences appear in all other probe brane cases 
as well. We will discuss these phenomena further later. 



4. Symmetry energy in quark matter system 

So far, we discuss symmetry energy in confined system. In this section we consider sym- 
metry energy in deconfined phase, namely, in quark matter instead of nuclear matter. The 
background geometry is a black brane. To introduce finite density or chemical potential, 
we introduce fundamental strings which connect black hole horizon and probe brane. The 
end points of fundamental strings attached on probe provide sources of U(l) gauge field. 
Since the fundamental strings can move freely along the D3 direction of the black brane 
horizon, the boundary system can be identified as the system of freely moving quarks. 

Similarly to the previous section, the symmetry energy can be understood as the energy 
cost to separate the number of up and down quarks from symmetric matter. From the D- 
brane point of view, we need to consider two probe branes, where different number of 
strings are attached so that the embedding of two branes are separated from each other. 
The schematic figure is drawn in Figure |5[ 

The boundary condition at the horizon is determined by the regularity of the equation 
of motion at the horizon; 

Y„ - = tan 6>, (4.1) 

where 6 is the angle of probe brane at the horizon. By substituting the embedding solution 
into (^1]), we can get symmetry energy as a function of density and quark mass. 
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(a) (b) 

Figure 5: Schematic figure of D-brane configuration for (a) symmetric matter and (b) asymmetric 
matter. 



4.1 D4 brane background 

To describe deconfined phase at the boundary theory, we introduce black hole geometry 
based on D4 brane; 



ds 2 = (|) 3/2 (f(u)dt* + dx 2 + dxl) + (^) 3/2 + u*dnl 

/ TT\ 3 / 4 O-irAT /7T„x 3 



2?riV r 

7l 



£4, /(C/) = l 



17/ 



R 3 = irg B N e ll (4.2) 



There is a horizon at J7 = Uq, and the Hawking temperature and the horizon radius are 

r-2 



related by [23 



9 



(4.3) 



Introducing a dimensionless coordinate £ by = n ^ (rn , the background geometry be- 



comes 



3/2 



3 /2 3 /2 2 

& 2 =(^) (/(f/Mt 2 + ^ 2 + ^) + (^ (|) (d£ 2 + £ 2 ^ 2 ), (4.4) 



where J7 and £ are related by 



3/2 



f/f) 



£3/2 



and / 



i-e 



-3\ 2 2 



(4.5) 



In this geometry, we rescale the coordinates such that the black hole horizon is located at 
£ = 1. The induced metric of probe Dq brane can be written as 

3/2 /,,2 



ds 2 Dq = (j^J £ 3/ V {^dt 2 + dx 2 }j 



R 3/2 



1/3 



£3/2 



(4.6) 



The form of free energy of D q brane can be obtained from ( A. ICQ as follows 



U). 



Q2+(£3/4^ + ) 



/ 2 V3\ n 
d-1 / p z u+ \ 



£3/4 



(4.7) 
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where n = q — d—1. By substituting the embedding solution into eq.( p.9[) , we can calculate 
symmetry energy in terms of density Y^, the asymptotic value of probe brane, which is 
related to quark mass and temperature by 



1< 

U 



9L 



m„ 



-m n 



S/T 



R 3 T 2 



(4.8) 



Therefore, if we fix quark mass, large value of Yqo corresponds to low temperature and 
small Yqo to high temperature. 

• D6 probe brane: 

In this case, at low temperature (large value of Y"oo); the symmetry energy S 2 = 0.5Q 1 / 2 , 
which is same as confining case. At high temperature (small value of l^o), symmetry en- 
ergy increase linearly in density. We can understand this from eq.(|2.9|). The symmetry 
energy is integration of a function from p m in to infinity. In the case of nuclear matter 
system, p m in is zero since the probe brane ends at the tip of baryon vertex. On the other 
hand, in deconfining phase, probe brane ends on black hole horizon with non-zero p m in- 
Therefore, embedding for small asymptotic value start near equator p m in = 1- If density is 
small(Q << 1), we can ignore Q 2 term in denominator and symmetry energy is just linear 
in Q. However, for low temperature, Yoo is small and p m in — > and therefore Q 2 in the 
denominator contributes to the integral by scaling the variables, giving the power behavior 
we observed above. The result is described in the figure |6](a). 

• D4 probe brane: 

In the case of -D4 probe brane, we get similar result in confining case which is lineally 
increase as density increase. See the figure 0(b). 





(a) 



(b) 



Figure 6: Density dependence of symmetry energy for several Y^. 



• D2 probe brane: 
For D2 brane case symmetry energy become 



So 



Q 



(Q 2 + 4)3/2 7 ^ 



dp 



2/3 



(l + Y 2 ), 



(4.9) 



which is also independent of details of embedding. And the result is almost same as Figure 

Hfc). 
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4.2 D3 brane background 

We now consider the black D3 brane geometry; 

U 2 / dll 2 \ 

ds2 = w ( f{u)dt2 + + r2 \WW + d v ' (4 ' 10) 

where R A = 2\l A s and f(U) = 1 - (U /U) 4 . Here, A = g\ M N c is the 't Hooft coupling of 
the YM theory. There is a horizon at U = Uq, and the Hawking temperature is given by 

T = -% = H2 . (4.11) 

Introducing a dimensionless coordinate £ defined by d£ 2 /£ 2 = dll 2 / (U 2 f (U)) , the bulk 
geometry becomes 

ds 2 = g (fdt 2 + dx 2 ) + ^ (de 2 + £ 2 dfi§) , (4.12) 

where £ 2 = y 2 + p 2 and /) is the radius of the 3-sphere. U and £ are related by U 2 /Uq = 
i(£ 2 + l/£ 2 ) and / = (1 — £ 4 ) 2 /(l +^ 4 ) 2 - The induced metric on Dq probe brane and 
Hamiltonian density can be written as follows 



'o; 2 



= T ? y— (l+y2)^Q2 +p 2n a; d. ( 4 . 13 ) 

With completely same method with the previous section, we calculate symmetry energy 
for each probe brane(L>7, D5 and .D3) in terms of density and asymptotic value of probe 
brane Y^. In this case, the relationship beetween the asymptotic value of probe brane and 
temperature(or quark mass) is given by 

Y -^L (4 14) 

The result is drawn in Figure [7[ Same as the previous section, we fix quark mass. Then, the 
asymptotic value of probe brane proportional to the inverse temperature. In the case 
of D7 brane, for low temperature(l^o = 15), the symmetry energy goes like S2 = 0.32Q 1 / 3 
while at high temperature, the symmetry energy grows linearly. The symmetry energy 
line at low temperature for D5 probe case goes like S2 = 0.5Q 1//2 and it becomes linear at 
high temperature. In the case of D3 probe, symmetry energy has linear behavior for all 
temperature. 

5. Scaling property and universality classes 

In this section, we want to understand the scaling behavior of the symmetry energy which 
is calculated in various models and various phases. We first tabulate all the results in table 
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Figure 7: Density dependence of symmetry energy for several Y c 



Table 1: DA brane background 
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Table 2: D3 brane background 
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1 and table 2. We emphasize that the result is based on the actual numerical calculation 
without any approximation. 

We now re-derive these power behavior of symmetry energy analytically. To do this, we 
consider the ideally simplified case: BPS metric background and flat embedding of probe 
branes. In this case, background geometry becomes geometry of black Dq brane; 

ds 2 w = Z- 1/2 (-dt 2 + dx 2 p ) + Zl l2 dx\, (5.1) 
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3-p 

with e 2<t> = Z P 2 .We now consider the general result for the symmetry energy 

S 2 = 2r q [ dp Qv^^^ , (5 . 2) 
Q2 + 4e -2^lG^ N ' 



which is already given in eq. ( |2.9[) . Then, the term in the square root becomes 1 and 

/ G±± \ q d 1 E=A _i q-d-l 

e 2 ^Ciz ( £2 ) = -^p 2 • Z p 2 • Z p 2 

= z| (p+9 " 2d - 4) . (5.3) 

The value of p + q — 2d is precisely equal to the number of Neuman-Dirichlet (ND) di- 
rection of Dp/Dq system. Therefore, if we focus on the system which is supersymmetric 
configurations or a smooth deformation of them, the exponent becomes zero. In this case 
the symmetry energy (|2.9| ) for flat embedding can be calculated analytically: 



Qp 



2n 



1 



Si = 2r q I dp- = c n Q™, (5.4) 



Q 2 + 4p 2 « 

where c n = 2r q r^/n 2 " anc ^ n = Q — d— lso that the density dependence of 

1 

symmetry energy is S2 ~ Qi- 11 - 1 . This result reproduces all the result we obtained in 
previous section numerically. 

Notice that both the background and embedding used here are far from the real sit- 
uation: real background is a deformation of such BPS solution and the embedding is 
non-trivial deformation from such a flat embedding. Nevertheless, the scaling exponent is 
the same as the actual configuration used for numerical computation of previous sections. 
The point is that neither smooth deformation of the metric nor the deformation of the em- 
bedding shape seem to change the scaling behavior of the symmetry energy. The exponent 
of symmetry energy depends only on the dimensionality of probe brane and dimension of 
non-compact directions. Therefore the scaling exponents depend only on the universality 
classes. 



6. Discussion 



In this paper, we calculated the asymmetry energy for both nuclear matter as well as the 
quark matter. The symmetry energy has a power like density dependence with charac- 
teristic exponent which is invariant under the smooth deformation of the metric as well 
as smooth deformation of the embedding. It only depends on the dimensionality of the 
D-brane system modeling the QCD dynamics. Therefore it is a index for the universality 
class. The physical interpretation of the scaling exponent is still open question but we give 
a trial interpretation below. 

Discussion: Non-fermi liquid nature of the nuclear matter. According to the fermi 
gas model of nuclei, the Energy of the nuclei is the sum of the kinetic energies of neutrons 
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and protons which are moving non-relativistically, namely, 

E = -(N p e Fp + N n e Fn ) (6.1) 

= ^e F A + ±e F A6 2 + 0(5% (6.2) 

where A = N p + N n is the mass number and e F is the fermi energy for A, 5 = |^ p + ^ n | . 
One can see that the symmetry energy per nucleon is e F /3. Therefore fermi gas model 
demonstrates the origin of the symmetry energy as the Pauli principle. What happens if we 
include the interaction energy of the nucleons? The answer is largely unknown. Depending 
on how one includes the interaction, the answers are different from one another. In some of 
the traditional approach, interactions are taken care of by adding a polynomial in density, 
which means that scaling property (~ p a ) should be destroyed by the interaction. However, 
this is not what we get. In our approach, scaling property remains even in the large 
density limit. Notice that for the most interesting case of D4/D6 and D3/D7, the scaling 
exponents are 1/2 and 1/3 respectively. Comparing this with the non-interacting gas which 
gives S2 ~ p 2//3 , one can see that the symmetry energy of DA/D6 implies an anomalous 
dispersion relations which is neither relativistic nor non-relativistic one. For the D3/D7 
case, the scaling is the same as the relativistic particles. If we naively extrapolate the 
relation S2 ~ e F , S% ~ p 1 ^ 2 means the anomalous dispersion relation e ~ A; 3 / 2 . One may 
want to associate this as the example of transformation of a particle to un-particle [^] 

One may also expect that such anomalous dispersion relation is related to the non- 
fermi liquid nature of the strongly interacting fermion system. In the presence of the 
fermi sea, one expect that elementary excitations are quasi particles with renormalized 
charge and mass. However, when the interactions are strong, such quasi-particles will lose 
applicability. Using the AdS/CFT duality and utilizing the AdS at UV and the AdS2 at 
the IR, it was shown in[^5| that 

G(u, k) = — 5-, (6.3) 

where k± is the momentum measured from the fermi surface. If v < 1/2, the dispersion 
relation becomes oj ~ k 1 [ 2h ' . For example, if v = 1/3 we get oj ~ k^[ 2 . In our case, the 
situation is more subtle since the fermi surface is not stable under turning on temperature 
and fermi sea is not so calm: the low lying states below 'fermi sea' can be all excited so that 
it should be called fermi ball |2(| rather than fermi surface, and k± should be replace by k. 
These phenomena are all beyond the fermi liquid behavior. For D3/D7 case, the exponent 
indicate that the system is marginally fermi liquid case. More systematic investigation on 
this matter is strongly desired. 
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A. Free energy of Dq brane 



We start from 10 dimensional metric for Dp brane with Lorenzian signature as 

dsj = -G tt dt 2 + G xx dx 2 + G rr dr 2 + G±±d£l 2 ± 

= -G tt dt 2 + G xx dx% + Gj_l (Jr^dr 2 + dQlj , (A.l) 

where dfl 2 ^ is metric of S , ( 8 ~p). To make the orthogonal space to Dp brane, we introduce 
new coordinate; 

§^dr = |, (A.2) 



then, the metric (A.l) becomes 



ds 2 = -Gudt 2 + G xx dx\ + ^ {de + fdfl 2 ± ) . (A.3) 

Now, we put probe Dq brane in this background with d dimensional uncompact direction. 
We can decompose transverse direction into parallel and perpendicular to the probe brane 
as follows 

ds lo = -Gttdt 2 + G xx {dx 2 d + dx 2 p _ d ) + -— - [dR 2 q _ d + dB%_ p _ q+d ) 



-Gudt 2 + G xx {dx 2 d + dx 2 

G -L-L ( iJ2 , Jjo2 , ,iv--2 : -1.-2. 



p-d) 



+^ (dp 2 + p 2 dtf q _ d _ l + dY 2 + Y 2 dn 8 - P - q+ d) ■ (A.4) 

We assume that only Y coordinate has p dependence and the other coordinates which 
perpendicular to probe brane is constant, then the induce metric on probe brane can be 
written as 

ds 2 Dq = -Gudt 2 + G xx dx 2 + ^(1 + Y 2 )dp 2 + < ^±p 2 dn 2 q _ d _ 1 

= -Gudt 2 + G xx dx\ + G pp dp 2 + GnndnJ-*-!, (A.5) 

where Y = dY/dp. 

To introduce number density in boundary theory, we turn on time component of U(\) 
gauge field At(p) on the probe brane and set all the other components to be zero. Then 
we can write DBI action as follows; 

S Dq = J da g+l e-^^det(g + 2vra'F) 

= -p q j dtdtxdnq-d-i dpe-ty/GixG^-^GttGpp - F 2 ) 



dtdpCoq, (A. 6) 

where F == 2-ira'F t p = 2-na' d p A t {p). 
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From the equation of motion for gauge field, we can define conserved charge; 
dC Dq _ T q e-<t>G d xx G^- l F _ Q 
dF ~ ■'G* a C%t\G tt G„-F*) W ' 



(A.7) 



where r q = p q VdO, q -d-i- Then we can get the relation between field strength and charge 
as 

F = QV* 7 ^ (A . 8) 

where Q = Q/2ira'T q . 

After Legengdre transformation for Lagrangian density, we can get free energy of D q 
brane; 

F Dq = J dp {^^j^ 1 ~ £o q 

= r q f dpy/G^^Q* + e-^GUG&t 1 (A.9) 

= r q J dp^G tt ^(l + ^Q 2 + e-^Gi x (^p^ 9 ' (A.IO) 
B. Baryon vertex and force balance condition 

In this work, we consider D3 or DA brane as a background. In D3 brane background, D5 
brane wrapping on S 5 with N c fundamental strings is interpreted as baryon vertex and 
a spherical DA brane on S 4 is baryon vertex in DA brane background. In the both case, 
background metric can be written in general form; 

ds\ Q = -G tt dt 2 + G xx dx 2 p + ^ (df + ednj_ p ) 

= -Gudt 2 + G xx dx 2 p + ^d£ 2 + G ±1 (d9 2 + sin 2 9dn 2 7 _ p ) , (B.l) 

here, we assume that the D{8 — p) brane has SO(7 — p) symmetry, then £ depends on polar 
angle 6 only. Then the induced metric on D(8 — p) brane is 

ds 2 BV = -Gudt 2 + + £' 2 ) d8 2 + G ±± sin 2 8dn 2 _ p . (B.2) 

The DBI action for D(8 — p) brane can be written as followes; 

Sbv = -fis-p J d 9 - p ae-^^det(g + 2ira'F) + p 8 _ p J 2<na' A t A G (7 _ p) , (B.3) 

where ps-p is a D-brane tension and G(7_ p ) is Ramon-Ramon field which couples to the 



original Dp brane. After substituting induce metric (B^), we can get 



Sbv = r 8 -p J dtdO sin( 7 -P) 6 -R^e^G~\ (£ 2 + £ /2 ) — F 2 + (7 — p)A t 
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= / dtdOCsv, 

where 

The displacement can be obtained by derivative the action with respect to F, 



.M, n R 7 -P, At = IvolAt. 



r-p 



dC BV T 8 - p RP- 7 sm 7 -P9e-< t> G ± 2 ± F 



dF 



GttGi 



and the equation of motion for gauge field is given by 

d e D(0) = -(7-p) sm 7 - p 6, 



-D{6), 



(B.4) 



(B.5) 



(B.6) 



(B.7) 



where D(9) = D(6)/t%- p . By integrating, we get the solution in terms of hyper geometric 
function as follows 



D(9) = c + (7 -p) cose 2 F X Q.E^, *,«*», 

For example, in the case of D4 brane (p = 4), we get 

D(0) = c + cos(9(3-cos 2 f9), 



(B. 



(B.9) 



which is same as result in 16]. The integration constant cq is determined such that D 
vanishes at 6 = which means all fundamental strings are attached on north pole of 
D(8 —p) brane. After substituting the solution of equation of motion, we can rewrite DBI 
action in terms of displacement which is called as 'Hamilotnian' because the procedure to 
get it is similar to Legendre transformation. 



Ubv = T~8- 



G tt G ±± ^ + ^/f)2 + R 2( P -7) Sin 2(7-P) Q e -2^ G 7- P> 



^ 2 



(B.10) 



As discussed in previous works, we can calculate force at the cusp of single D(8 —p) brane 
as follows, 



Fbv 



sn 



BV 



dU 



BV 



on— shell 



(7- P )v^r(V 



2- 1 :[5=2 



GuG 



N c G U G±± i' c 



(B.ll) 



where £ c and £' c denote to the value of £ and it's derivative at the cusp. The overall factor 
becomes N c /2ira' for both of D3 and D4 case. 
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Now, we can get boundary condition for probe Dq brane by imposing 'force balance 
condition'. The force at the cusp of probe brane can be obtained from Hamiltonian of 
probe brane ( |A.1Q ); 



F, 



Dq 



Dq 



dU 



Dq 



on— shell 



dY 



p=0 



Q G tt G±± Y c 



2na' 



l + Y- 



(B.12) 



where Y c and Y c denote to the value of Y and it's derivative at the cusp of probe brane. 
From the force balance condition 



F Dq (Q) = N B ■ F BV = ^-F BV 
we can get the boundary condition for probe brane 

V — 2£. 

C ~Y 



(B.13) 



(B.14) 



which does not depends on type of original brane or probe brane. 
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